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(GRBs) in the BATSE 4B catalog can be used to constrain the 
average redshift < z > of the GRBs. In this paper we investigate 
p3 ■ the effect of evolving lenses on the < z > of GRBs in different 

cosmological models of universe. The cosmological parameters 
Q and A have an effect on the < 2 > of GRBs. The other 
factor which can change the < 2; > is the evolution of galaxies. 
We consider three evolutionary model of galaxies. In particular, 
we find that the upper limit on < 2; > of GRBs is higher in 
evolving model of galaxies as compared to non-evolving models 
of galaxies. 
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1 Introduction 

The use of gravitational lensing as a tool for the determina- 
tion of cosmological parameters (e.g. Hq, Qq, Aq) has frequently 
been discussed. M To constrain these parameters either QSOs or 
galaxies have been used as sources. To use Gamma Ray Bursts 
(GRBs) as a source for gravitational lensing is not a new idea. 
As first discussed by Paczynski,^'!^ if GRBs are cosmological then 
they should be gravitationally lensed just as quasars. 

There is now overwhelming evidence that the majority of 
GRB sources lie at cosmological distances. The detection of 
GRB990123 which is believed to lie between redshift^ 1.6 < z < 
2.14, the identification of the host galaxyi for GRB981214 at z = 
3.42 and the detection of afterglows of many others GRBs^'M in 
X-rays, optical and radio domain by BeppoSAX satellite are all 
suggestive of GRBs being at cosmological distances. Therefore 
it is expected that the GRBs must be gravitationally lensed as 
the probability of lensing increases with the source redshift. The 
primary effect of a gravitational lens on GRB would be to create 
more than one image of the burst. These images could not be 
angularly resolved with the present technology since BATSE's 
spatial resolution is insufficient to resolve the images angularly 
but they could be temporally resolved. Several authorsli^l^^ have 
presented detailed calculations of GRB lensing which can be 
used as a probe in the search for dark matter in the form of 
compact objects. 

So far no lensed GRBs have been detected with BATSE which 
clearly implies an upper limit to the average redshift < z > of 



GRBs as first calculated by Holz, Miller & QuashnockUl [HMQ]. 
HMQ calculated an upper limit to < 2; > which is independent 
of the physical model of GRBs. They have calculated < z > for 
GRBs in different cosmological models. 

Till date a substantial fraction of all GRBs observed occur 
at significant cosmological distances ( z ~ 0.8 — 3.4). Further 
the " no host galaxy problem" pushes the GRBs to be, either at 
very high redshift { z > 6) or not to be in normal host galaxies.li^ 
Several authorgi3'EJ propose that the GRBs rate should trace the 
star formation rate in the universe and consequently it places 
the very dim burstHi at z > 6. HMQ calculated the < z > of 
GRBs in a non-evolving model of galaxies (lenses) and the < z > 
at 95 % confidence level (CL) comes out to be < 2.2, 2.8, 4.3 for 
different combinations of (Q, A) namely (0.3, 0.7), (0.5, 0.5) and 
(0.5, 0.0) respectively. These results get considerably modified 
when evolution in the properties of lensing galaxies (lenses) is 
considered. The modified results permit a larger value of < z > 
and are thus consistent to some extent with the high - z GRB 
scenarios. The purpose of this paper is to put the upper limit on 
< z > of GRBs in different models of evolving lensing galaxies. 

We refine the analysis of HMQ by considering the evolution 
of galaxies. Normally the comoving number density of galaxies 
is assumed to be constant while calculating lensing probability. 
But it is an oversimplification to assume that galaxies are formed 
at a single epoch. Evolution tells us how the mass or number 
density of the lens varies with cosmic time scales. Merging be- 
tween galaxies and the infall of surrounding mass into galaxies 



are two possible processes that can change the comoving density 
of galaxies and/or their mass. The effect of galaxy merging or 
evolution has been studied by many authors. li^T^S Most of them 
have focused on the statistical properties of gravitational lenses 
and the limits on the cosmological constant. This work is an 
attempt to check how galaxy evolution changes the upper limit 
on the average redshift of GRBs. 

The structure of the paper is as follows. In section 2 we 
describe the evolutionary model of galaxies that we use. The 
lensing probability with evolving lenses is given in section 3. 
In section 4 we put the upper limits on the average redshift 
of GRBs, assuming that no lensing events are present in the 
BATSE 4B catalog. Results on < z > of GRBs with evolving 
lensing galaxies in different cosmological models are described 
in section 5. Discussion is given in section 6. 

2 Evolution Of Galaxies 

Galaxy mergers and gravitational lensing are interlinked with 
each other as the merging of small galaxies gives rise to elliptical 
galaxies which further act as gravitational lenses. We consider 
three evolutionary models of galaxies 

Fast Merging 

This evolutionary model was proposed by Broadhurst, Ellis 
& Glazebrook (1992) [BEG]E3 in order to resolve the faint galaxy 
population counts. The basic point of this model is to introduce 
strong number evolution or to assume that galaxy numbers are 



not conserved with cosmic time but the total comoving mass 
remains conserved. This model assumes the comoving number 
density of the lenses to vary with cosmic time as : 

n{6t) = f{6t)no (1) 

where 5t is the look-back time and the subscript '0' indicates 
the present day values. The characteristic mass ( or disper- 
sion velocity ) for self similar galaxy mass function of Singular 
Isothermal Sphere (SIS) lenses at 6t is 

v{6t) = [f{St)]-\ (2) 

This model is better than the model proposed by Volmerange 
and GuiderdonicB where the merging goes exponentially with 
(1 + z), because the merging rate doesn't become high at early 
time as explained by BEG. According to the BEG model if we 
had n galaxies at look-back time 6t each with velocity dispersion 
■u, they would by today have merged into one galaxy with a 
velocity dispersion [f{6t)]v. The function f{6t) describes the 
time dependence and the strength of merging : 

f{St) = exp{QHo6t) (3) 

where Hq is the Hubble constant at the present epoch and Q 
represents the merging rate. We take Q = 4 as suggested by 
BEG.I^ The look back time 6t is related to zi through 



where F{y) = Qo(l + vf + (1 - Oq - Ao)(l + y^ + Ao 

where Qq = SttGpo/SHq , po is density of matter, Aq = 87TGpyo/3HQ 

and pyQ is density of vacuum 



Slow Merging 

In this slow merging modell^ the mass of an individual galaxy 
increases with cosmic time as t'^'^ while the comoving number 
density varies with cosmic time as t^'^'^, so the total mass of 
the galaxies within a given comoving volume is conserved The 
cosmic time t starts from the big bang. We also assume that the 
power law relation between the mass and velocity is M oc v^'^ 
for elliptical galaxies.^ Then 



n{5t) = TiQ 



1 



6t 
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v{St) = vo 
where to is present age of the universe 






(5) 



(6) 



Mass Accretion 



Mass accretion is the key factor for evolution of galaxies. In this 
model^ the comoving number density of the galaxies remains 
constant but the mass of galaxy increases with cosmic time as 
t^'^ in the same way as in the slow merging model. 



n{5t) = no{constant) (7) 



v{6t) = vq 
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(8) 



3 Basic Equations For Gravitational Lensing 
Statistics 

The differential probability dr of a beam encountering a lens in 
traversing the path of dzi is given byc^'B^ 

cdt 
dr = nL{z)a--—dzL, (9) 

dZL 

where nL{z) is the comoving number density. 

The Singular Isothermal Sphere (SIS) provides us with a rea- 
sonable approximation to account for the lensing properties of a 
real galaxy. The lens model is characterized by the one dimen- 
sional velocity dispersion v. The deflection angle for all impact 
parameters is given by o; = ^t^v^ jf? . This lens produces two 
images if the angular position of the source is less than the crit- 
ical angle i5cr^ which is the deflection of a beam passing at any 
radius through a SIS: 

4, = aD^sl^os. (10) 

Here we use the notation Dql = d{0, z^), D^s = d{zL^ zs), Dqs = 
(i(0, zs), where d{zi, Z2) is the angular diameter distance between 



the redshift zi and Z2-^ Then the critical impact parameter is 
defined by acr = DoiPcr and the cross- section is given by 



The Evolutionary Model 

The differential probability dr of a lensing event in an evolu- 
tionary model can be written as: 

leTT^ , 4-r. / 4 

—-^(j)^v^r[a + - 
cH§ V 7 

n , ^i(DoLDLs\^ 1 cdt 



Av{^.^l^A fa^^(l-|) 



dT = —-^(j),vtr(a + - + i]f{6ty 






where F = ^^ 0^ i;^r ( a + - + 1 ). By substituting the values 

'' V 7 pZ_ 

of 0*, a as given by Corray et al.^ and "U*, 7 as given by Naka- 
mura & Suto^ we get F = 0.035. where f{6t) = exp{QHQ6t) 

for fast merging and /((5t) = [1 — — ] for slow merging. In 

case of mass accretion /(^t) = f 1 — ^j but the exponent of 
/((5t) for mass accretion model in eq. (12) becomes (— 1 — -) as 
the total mass in galaxies increases with time. 

The Non Evolutionary Model 



In the non merging model the optical depth is given byl^ 



(DolDls\^ 1 cdt 
4 Bound On Average Redshift Of GRBs 



Adopting the simplest matter distribution for the lensing galaxy 
as a Singular Isothermal Sphere (SIS), we calculate the proba- 
bility t{z) of a beam from a source at redshift z is imaged by a 
lens in the filled beam approximation. Apart from the cosmo- 
logical parameters H and A, the lensing probability depends not 
only on the model used for describing the evolution of galaxies 
but most importantly on the parameter F.E^'Ey'E3'l^ F parame- 
terizes the distribution of galaxies as well as their effectiveness 
in the lensing process (defined in section 3). Thus the elliptical 
galaxies contribute strongly while spirals have a negligible ef- 
fect. The Schechter distribution function often forms the basis 
of these estimates. However it is being increasingly realised that 
sub-distribution of specific types of galaxies have to be further 
taken into account. We take F = 0.035 which lies in the current 
estimate rangeHiUl'il 0.02 < F < 0.05. In the SIS model, 
the lensing event always consists of two images and the third 
central image is too faint to be observed. We take a constant 
BATSE efhciencyi^ e = 0.48 to see either image and e^ = 0.23 to 



see both images. So far no multiple images has been detected, 
which may due to the low instrumental efficiency. 

In order to put the bound on the < z > of GRBs, we fol- 
low the methodology of HMQ. As the distribution of GRBs in 
redshift is unknown, we use the simplifying assumption that all 
the sources are at the same redshift. We estimate the number 
of image pairs as 

iV<z> = {Ntot/e).e\r{< z >) = iVto,er(< z >) (14) 

where Ntot be the total number of observed bursts in the 
BATSE 4B catalog (which are 1802) then Ntot/^ are actual burst 
sources above the BATSE threshold. The < z > m. the above 
equation tells the average redshift of GRBs, assuming that all 
the sources are at this redshift. If we define ^^m as the bright- 
ness threshold below which identification of lensed images from 
light-curve comparison is impossible, then the expected number 
gets modified tot^l^ 



N. 



<z> 



roo 1 

= N,oteT{< z >) J dB{<p) (15) 

•^0 [1 + (pmin/(p{z)r 



We plot this N^z> with < z > in different cosmological mod- 
els as shown in Fig. 1, Fig. 2 and Fig.3 respectively. Here B{(j)) 
is the the observed BATSE brightness distribution and the in- 
tegrand is the conditional probability that the both images are 
above the brightness threshold. The value of this integral is 
equal to 0.57.0 
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5 Results 

We have taken three representative values of {Q, A) as (0.3, 0.7), 
(0.5, 0.5) and (0.2, 0). With these values we calculate the ex- 
pected number of observable image pairs in the BATSE 4B as a 
function of average redshift for the mass accretion, fast merging, 
slow merging and no evolutionary models. We find that the best 
limit arises with a large cosmological constant where the lensing 
rate is quite high. At the 95 % confidence level, we find an upper 
limit on {z) < 4.3, 7.8, 9.8 in the non-evolutionary model for 
(n. A) values of (0.3, 0.7), (0.5, 0.5) or (0.2, 0) respectively. The 
upper limit on {z) at 95% CL in other evolutionary models are 
described in table 1. At 68 % confidence level, the slow merg- 
ing model gives upper hmit on (z) < 2.3, 3.3 or 3.9 for (n,A) 
values of (0.3, 0.7), (0.5, 0.5) or (0.2, 0) respectively. Similarly 
fast merging model gives (z) < 2.9, 4.3 or 5.2 for (Q, A) values 
of (0.3, 0.7), (0.5, 0.5) or (0.2, ) respectively as shown in fig.l, 
fig. 2, and fig. 3 respectively. The mass accretion model at 68% 
CL gives an upper limit on (z) < 6.1 only for {Q = 0.3, A = 0.7). 
The result is described in Table 1 and Table 2. 

6 Discussion 

We notice that the evolution of lensing galaxies serves to increase 
the < z > of GRBs as compare to non- evolutionary model of 
galaxies. 

Throughout the paper, we have used the filled beam approxima- 
tion (standard distance)^^ in which "smoothness parameter" 
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a ( which measures the degree of inhomogeneity of the universe) 
is equal to 1. On the other hand with a = one gets empty 
beam distance in which there is negligible intergalactic matter 
and the hue of sight to a distant object doesn't pass close to 
intervening galaxies. 1^'^ At present we don't know what value 
of smoothness parameter in the distance formula describes our 
universe best. Ehlers and Schneiderc^(1986) argued that the di- 
rection to the source cannot be a random variable in the clumpy 
universe. At lower redshift the probability of lensing in different 
distance formulations remains approximately the same while at 
higher redshift the lensing rate become different in different dis- 
tance approximations. So the question arises, which is the best 
distance approximation? But recently in an elegant paper, Holz 
& Wald (1998)cll have exhibited a new method for calculating 
gravitational lensing rates. This method is free from the ambi- 
guities present in the distance formula. In this formulation, the 
optical depth at any redshift is greater than the optical depth 
calculated using the angular diameter distance in the filled beam 
approximation. Thus we expect that this formulation will in- 
crease the expected lensing rate, and hence decrease the upper 
limit on < z > of GRBs. 

Another reason why our lensing rate could be an underesti- 
mate because of the assumption that lensing is only due to the 
SIS galaxies. But there are other lens models also like: point 
like mass as deflector, isothermal sphere with a softened core, 
asymmetric lens , spherical model (e.g. King model), cosmic 
strings and constant density sheet which in principle could also 



12 



contribute. Many authorsll3c5E^ have explored the possibihty 
of GRBs lensing by point masses. Grossman and Nowak^ also 
studied the lensing of GRBs by asymmetric and non singular 
isothermal spheres lenses. These additional lensing objects may 
change < z >. 

The overall probability of strong lensing depends directly on the 
parameter F, which further depends upon four parameters a, 7, 
-Uh., 0*. There are many uncertainties associated with these pa- 
rameters as mentioned in sec. 4 and discussed by several au- 
thors.l^'I^E§ED These uncertainties can change the value of F 
by upto 30% and hence affect our results. HMQ used the value 
of F = 0.1 which is nearly three times higher than the value of 
F used in this paper. Therefore the upper limit on the < 2; > of 
GRBs in HMQ paper is less than our value of < 2 >. 

The ability of BeppoSax satellite to detect GRBs opens up 
a new era in the studies of GRBs. Many counterparts of GRBs 
have been detected in X-ray, optical and radio domains from 
which important insight in this field has been gained. Never- 
theless, the nature of the central engine that accelerates the 
relativistic flow is still not very clear. In order to calculate the 
amount of total energy in the burst, the precise measurement 
of redshift is required. For example, spectroscopic observations 
show that the host galaxy for GRB971214 is at a redshift z = 
3.418. Given this high redshift and the known fluence of this 
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GRB the 7 - ray energy release of this burst is unexpectedly 
large, about 3 x lO^^erg, assuming isotropic emission.^ Energy 
released in other forms of radiation is not included in this en- 
ergy calculation. The currently favored model for GRBs is coa- 
lescence of neutron stars.ea The coalescence model is expected 
to released about lO^^erg in the form of electromagnetic energy. 
So there is big gap between the model based calculation and the 
parameters based calculation. Similarly, with the combination 
of a redshift > 1.6 and fluence for GRB990123i^ imply that 7 
ray energy release in this burst is 3.4 x lO^^erg, assuming the 
emission is isotropic. This would again put strain on any GRB 
model based on merging of neutrons star or black holes. At 
higher redshift this problem becomes more severe. We may be 
forced to consider even more energetic possibilitied^ or to find 
ways of extracting more electromagnetic energy in coalescence 
models or we have to search for ways to reduce the estimated 
energy^ ^ release from GRBs to resolve this problem. 
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Figure Captions 

Figure 1. Expected number of observable image pairs in the 
BATSE 4B {N^z>) as function of < z >, the effective average 
redshift of GRBs, in fiat universe with Q = 0.3 and A = 0.7 . 

Figure 2. Expected number of observable image pairs in the 
BATSE 4B {N^z>) as function of < z >, the effective average 
redshift of GRBs, in flat universe with il = 0.5 and A = 0.5 . 

Figure 3. Expected number of observable image pairs in the 
BATSE 4B {N^z>) as function of < z >, the effective average 
redshift of GRBs, in an open universe with Q = 0.2 and A = 0.0 
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Table 1. Limits On The Average Redshift {z) Of GRBs In Different Models 
Of Galaxy Evolution At 95% Confidence Level with 1802 bursts 



VL + ^ 


No Evo. 


Slow Merging Fast Merging 


Mass Accretion 


0.3 + 0.7 


4.3 


5.0 7.6 





0.5 + 0.5 


7.8 


10.0 





0.2 + 0.0 


9.8 









Table 2. Limits On The Average Redshift (z) Of GRBs In Different Models 
Of Galaxy Evolution At 68% Confidence Level with 1802 bursts 



Q + A 


No Evo. 


Slow 


Merging 


Fast Merging Mass Accretion 


0.3 + 0.7 


2.2 


2.3 




2.9 6.1 


0.5 + 0.5 


3.0 


3.3 




4.3 


0.2 + 0.0 


3.5 


3.9 




5.2 
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Q = 0.3, A = 0.7 




Mass Accretion 
Fast Merging 
Slow Merging 
No Evo. 



Figure 1: 
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^ = 0.5 , A= 0.5 




Mass Accretion 
Fast Merging 
Slow Merging 
No Evo. 



95 % CL 



< z > 



Figure 2: 
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Q = 0.2, A = 



2 - 



Mass Accretion 
Fast Merging 
Slow Merging 
No Evo. 




< z > 



Figure 3: 
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